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We analyze the collective spin noise in interacting spin systems. General expressions are derived 
for the short time behaviour of spin systems with general spin-spin interactions, and we suggest 
optimum experimental conditions for the detection of spin squeezing. For Ising models with site 
dependent nearest neighbour interactions general expressions are presented for the spin squeezing 
parameter for all times. The reduction of collective spin noise can be used to verify the entangling 
' ' ' powers of quantum computer architectures based on interacting spins. 
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^ . I. INTRODUCTION 

' Parallel with the study of non-classical squeezed states of electromagnetic radiation increasing attention has 
I been devoted to the study of atomic spin squeezed states [p|-[l4[. In Ramsey spectroscopy a sample of N two- level 

T-H ' atoms is represented by the collective operators = X^iLi {oi — x, y, or z), where ja,i — aa,i/2 and aa_i are the 
^ Pauli operators for the z*'' atom. Atomic spin squeezed states are quantum correlated states with reduced fluctuations 
in one of the collective spin components, and they have possible applications in atomic interferometers and high 
\ precision atomic clocks. Wineland et al. [g| have shown that the frequency resolution in spectroscopy depends on the 
' spin squeezing parameter 
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f*) i] where _L denotes a direction perpendicular to the mean spin. The inequality < 1 indicates that the system is spin 
I ' . i: squeezed, and it has been shown that any state with < 1 is an entangled state p^JTl[ |. 

In 1993 Kitagawa and Ueda showed that spin squeezing is produced by simple nonlinear spin Hamiltonians in 
^ , analogy with the non- linear Hamiltonians leading to squeezed light |^]. A number of experimental proposals for 
atomic spin squeezing have appeared involving interaction of atoms with squeezed light ||^-|^ , quantum non-demolition 
measurement of atomic spin states and atomic coUisional interactions P,p^, and recently the first experimental 
realizations of spin squeezing have been achieved 

There is a link between the theory of quantum computing and the theory of spin squeezing through the observation 
that the register of quantum bits in a quantum computer constitutes an ensemble of effective spins, and a general 
' purpose quantum computer can obviously produce a spin squeezed state of its register qubits. Moreover, proposals 
for quantum computing may be partly implemented, e.g., with only a restricted class of operations, which do not 
suffice for general computing purposes, but which do lead to massive entanglement and possibly to spin squeezing. 
We have previously considered such reduced instruction set (RISQ) quantum computers with ions and atoms and 
we have pointed out that apart from their ability to address particular physics problem such as anti-ferromagnetism, 
they can also be used to synthesize useful quantum states. In this paper we present a theory for the spin squeezing 
expected for different models of interacting spins. These models encompass a number of theoretical proposals for 
quantum computing, and although spin squeezing may not be a particularly relevant property in, e.g., a spintronic 
or a quantum dot realization of a quantum computer, we wish to point out the possibility of verifying the entangling 
powers of the permanent or controllable interactions in these systems by simple measurements on the entire system. 
The detection of spin squeezing may constitute a useful diagnostic tool in the early stages of the construction of a 
quantum computer. 

The treatment of the most general case of interacting spin systems is a formidable task, dealt with by a number of 
ingenious approximations in the theory of magnetism in solid state physics. In section |^ of the paper, we consider 
the case of an initially known state of the spins, subject to an arbitrary interaction Hamiltonian. By application of 
Ehrenfest's theorem, we can determine the short time behaviour of the system analytically, we can identify the spin 



squeezing signal, and we can devise the optimum conditions for detecting this signal. In Section [II, we consider a 
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general Ising type Hamiltonian with the spins constituting a chain with only nearest neighbour interactions of varying 
magnitude. We obtain an analytical expression for the spin squeezing of the system at any future time, and we provide 
different examples of spin chains with constant, alternating and random couplings. 



II. SHORT TIME BEHAVIOUR FOR GENERAL PAIRWISE INTERACTIONS 

Consider a collection of spin 1/2 particles which only interacts through pairwise interactions. The most general 
Hamiltonian describing this situation is given by 

H^Y.^,^-m''-l^ (2) 

where denotes the transpose jk = {jx,k, jy,k, jz.k), and where rri'^' are real 3 by 3 matrices. We assume that the 
particles are initially prepared in a state where all spins are pointing in a given direction prepared for instance by 
optical pumping techniques, and we compute the time evolution of the noise in a component perpendicular to the 
direction of the mean spin. A convenient representation of this situation is by a collective spin state R{a, (3, 7)! J J), 
where R\s a, rotation operator given by the three Euler angles a, (3, and 7, and where \ JM) denotes an eigenstatc of 
J ^ and Jz with eigenvalues J( J + 1) and M . The noise in a direction perpendicular to the mean spin can be found 
by calculating the square of the rotated operator Jj^ — RJ^w . 

With the general Hamiltonian (|^) it is not possible to calculate the full time evolution of the noise. It is however 
straightforward by Ehrenfest's theorem to determine the short time evolution 

^^{AJ^r = ^{JJ\[H,Ji]\JJ), (3) 

where we have introduced the transformed Hamiltonian H = W HR and we have used ( J_l)=0. By using wJR = R J, 
where R is the 3 by 3 matrix representation of the rotation in coordinate space, we may express the transformed 

Hamiltonian as H = 'J2k^i jk^ ' ' i'' where m''' = R^ ■ m*^' • R. By calculating the commutator and evaluating the 
expectation value we find 

|(AJ.)^ = ^(0,1,0).M-| |. (4) 

In this expression we have introduced a matrix M_ = '^k^iUl'^''' , and transformed it with the rotation matrix M = 

In Eq. (||) the terms involving the kih and lih particles are jk^ •m'^' ■ ji and ji^ -m}^ ■ jk. The sum of these terms can 

also be written as l/2(jfe^ • (to*^' + m^'' ■ ji + ji^ ■ {rn}'' +rn'''- '"") ■ jk)- From this expression we see that without loss 
of generality we can assume the matrix M_ to be symmetric and hence diagonalizable and in a convenient coordinate 
system we have 



M=\ My \. (5) 

Mz 

By calculating M_ using the form (|^) for the matrix M and the well known rotation matrices R, we can find the 
time derivative of the noise perpendicular to the mean spin for any orientation of the spin by using Eq. (^. Since we 
start out in a state with = 1 and since d/dt{J )^ = at t = 0, the interaction produces spin squeezing if we can 
find any set of Euler angles a, /?, and 7 which gives a negative derivative in Eq. (jj). The optimal orientation of the 
spin is found by minimizing the derivative with respect to the angles a, (3, and 7. We find that the extrema of 
Eq. (^) are always with the mean spin along one of the eigenvectors of M. If the spin is polarized along the z-axis 

the change in the noise is maximal for the perpendicular components J±7r/4 = ^l"/^{Jx i Jy) and we find 



|(AJ±,/4)2 = ±^(M,-A4). (6) 
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Similar expressions are found if the spin is oriented along the x or y-axes. 

Note that no assumptions about the values of the coupling matrices are made, they may vary randomly for any pair 
(fc, I) of spins. Only the sum needs to be specified to determine the short time spin squeezing. The above argument 
only states that some squeezing will be produced, it does not say anything about the maximum squeezing. In the 
following section we shall analyze special cases of the interaction (||) where the squeezing can be calculated exactly 
for all times. 



III. ISING CHAIN WITH NEAREST NEIGHBOUR COUPLING 



Consider now a general model with arbitrary coupling constants between nearest neighbours in an Ising spin chain. 
The chain consists of N spins with the Hamiltonian 

N 

H = fi^Xtjx,tjx,i+i, (7) 

i=l 

where we identify the N + 1** spin with the first one in the chain. Depending on the value of xn, the chain can be 
equipped with open or closed boundary conditions. 

This Hamiltonian arises in recent proposals for quantum computation with atoms in optical lattices [ p^ , p^ . In 
these proposals the atoms interacts with the nearest neighbours and it has been shown that the interaction can be put 
in the form (Q) and that this interaction produces spin squeezing . Other application of this Hamiltonian in optical 
lattices can be found in |p8| . Spin squeezing is much less experimentally challenging to produce and to verify than 
the full quantum computer, and squeezing can provide a demonstration of the entangling capabilities of the setup as 
well as having practical application in atomic clocks. From the discussion in section |l| follows that for short times 
the optimal squeezing is produced by having the spin initially polarized along z and by looking at the noise in one 
of the components J±^/4 — l/V^{Jx ± Jy), where = TjJ2iXi- For longer time intervals the optimal 

noise reduction is in a component Jg = cos{9)Jx + sm{9)Jy with 9 ^ ±7r/4 but for simplicity we shall only consider 
9 = ±7r/4. 

The unitary time evolution operator for the Hamiltonian H can be written 

N 

U(t) = e-'"''^' = Jl ^^^tj.,^3.,^ + l ^ (g) 

i=l 

and the operators ii{t) in the Heisenberg picture are readily obtained 

/I \ f jxM 

jy,i{t) ] [ ^ COs{Xi~ltjx,i^l + Xitjx,i+l) -Sm{xi-ltjx,i-l + Xitjx,i+l) ■ jy,i(0) I . (9) 
jz.i^t) ) \0 sin(Xi-lija;,i-l +Xitja;,-i+l) COS (Xt-ltjx.i-l + Xitix,t+l) ) \jz.i{^) 

We assume an initial state where all spins are pointing up, i.e., the jz^i = \ eigenstate for all i, and Eq. (^) then 
provides the mean values of the collective spin components at later times: 

{Jx) = {Jy) = 

^j^)-\y.-o44)cos{^]. (10) 

Since the collective operator Jx commutes with H, it follows that ( Jj) retains its original value at i = 0: (Jj) = ^. 



2 ^ \ 2 y \ 2 



In order to calculate {Jy), we need to know the the expectation values {jy^i{t)jy^i+k{t)) ■ We always have {jy,i{t)'^) — \, 
and by direct calculation we see that {jy^i{t)jy^i^i{t)) = 0., whereas for fc = 2 we obtain 



{jyAt)jy,^+2it)) = \ cos (^^^y-^ sin sin (^^y-^ cos (^^^^ ■ (H) 

From Eq. (^) we see that {jy,i{t)jy^i+k{t) vanishes for k >3, and we thus obtain 

/ 72\ ^ , 1 f Xi-it\ ■ f Xtt\ ■ (Xt+it\ f Xt+2t\ 
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and similarly we can determine the expectation value 



1 ^ 

(JyJx + JxJy) ^ y ] si 



i=l 



Xi + Xi+1 



(13) 



The degree of spin squeezing can be characterized by the squeezing parameter for a particular direction of the 
noise reduction. With the mean spin in the z direction, and considering the reduction of the noise in the component 
Jg the squeezing parameter can be calculated by 



i 



2 _ COS2(0)( J2) + Sin2(0)(j2) + sin(0) COs{e){JxJy + JyJx) 



(14) 



and for the particular component Je=7!/A the resulting squeezing factor is 

+NY.ti cos (^) sin (^) sin (^) cos (^) - NY.ti sin (^^^f^i) 



C7r/4 — 



Ef=iCos(M)eos(2-±il) 



(15) 



which is a general expression for any set of values for the coupling coefficients Xi- For & uniform closed chain (all Xi 
identical), the above equation reduces to 



2 _ 1 + 0.25sin2(xt) - sin(xi) 
" cos4(xV2) 



(16) 



as found in dj. 



A. Ising model with few spins 



In the above expressions, coefficients Xi+fc with i + k > N assume in a cyclic manner the values of the coupling 
among the first spins in the chain. It should accordingly be noted that Eq. (^5|) is not valid if < 4 , for which special 
expressions apply. 

In the Ising model with only two spins H2 = 2fi,xjx,ijx,2, the related expectation values are obtained as (J^) — 
(Jy) — 0, {Jz) — cos (xt) , (J2) = (Jy) = 1/2, and {JxJy + JyJx) = — sin {xt) , from which we obtain the squeezing 
parameter 



sin (xt) 



cos^ (xt) 

The system can be squeezed, and the maximum squeezing ^^^^ = 0.5 occurs when t = it/ (2x) 
For the Ising model with three spins = ?iEi=i Xijx.i3x,i+i: the squeezing parameter is 

9 + 3 Ell sin (^) sin (^) - 3 Eli sin 



(17) 



" COS (M) cos 

and for a uniform closed chain this equation reduces to 
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(18) 



l + sin^ (f ) - sin {xt) 



/4 



3S4 (f ) 



(19) 



which is different from Eq.(p^). 



Finally, in the Ising model with four spins Hn — ftEj=i Xiix,ijx.i+i- The squeezing parameter is given by 

K(^J> + ^'»-4EtiSin(^^t) 



/4 



Ell COS (^) COS (2^) 



(20) 
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where 
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(J!) + (4) - 2 + i E sin (f ) sin (^) cos (^) cos ) . (21) 



For a uniform chain Eq. (|2C|) actually reduces to Eq.(16), but for general coupling constants Xi, Eq. ( |15| ) is only valid 
for TV > 5. 

Fig. 1 is a plot of the squeezing parameter as a function of time for few spins in a uniform chain. It shows that 
maximum squeezing occurs for N — 2 and that the degree and the temporal range of squeezing are both small for 
= 3 in comparison with the other cases. 
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FIG. 1. Squeezing parameter as function of time in the Ising model with uniform coupling coefficients. The solid curve is for 
N = 2 spins, the dashed cuve is for TV = 3, and the dotted curve is for A'' = 4 spins or more. 



B. Dimerized and random chains 



It is interesting to study a dimerized Ising chain with even number of spins N = 2M, where the coupling constants 
are chosen as Xi — x[l + (— l)*^^'^]- We have two values of the couplings Xo ~ + ^) for odd i and Xe — " ^) 
for even i. Such systems appear naturally in heterogeneous structures where every second site is occupied with one 
type of qubit/particle, and their mutual communication is provided through intermediate particles acting as short 
range data-bus elements in e.g. the quantum computer. Two overlapping optical lattices with two different kinds of 
atoms or atoms in two different ground states may be moved relative to each other in order to establish the dimerized 
chain Hamiltonian. 

For the dimerized Ising chain, Eq.(|l5|) reduces to 

2 _ 1 + 0.25 sin (xot) sin (xe^) - sin {xt) 
W4- cos2 (2^) cos2 (2^) ' 



and we see that the spin squeezing does not depend on the number of atoms. Obviously Eq.(|22|) reduces to (16) in the 
limit (5^0, and in the limit (5 — > 1 it gives Eq. (^^. Fig. 2 is a plot of the squeezing parameter as a function of the 
parameter 6. We see that both the squeezing range and degree of squeezing increases as 6 approaches unity, i.e., the 
introduction of a dimerized coupling makes the squeezing better. For S larger than unity, the denominator vanishes 
within the time interval displayed in the figure, as shown for S = 1.1 which causes a divergence of the squeezing 
parameter at = 7''/2.1 = 1.496. 
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FIG. 2. Squeezing parameter as function of time in the Ising model with dimerized coupling coefficients. The curves are 
obtained with the following values of S: (solid curve) , 0.5 (dashed), 0.75 (short-dashed), 1 (dotted), and 1.1 (dot-dashed). 
The curves with S — and S = 1 are recovered in Fig. 1. 



Let us finally study a random chain model in which spins interact with their neighbours with a fixed coupling 
constant but only with probability p, i.e., any coupling constant is x with probability p and zero with probability 
1—p. This model corresponds, e.g., to the atomic lattice system with a non-unit filling fraction, so that empty lattice 
sites appear. Different assumptions may be made about the correlations between different lattice sites, but let us 
for simplicity consider the case where the occupancies are completely uncorrelated. Simulations for such a model 
were presented in Ref. but we can now present analytical results (for simplicity presented only for the 9 — tt/A 
direction). 

We introduce /i = in terms of which we can write the mean values of trigonometric functions. 



cos( — ) 
sm( — ) 



pcos + (1 — p) 



p sm n 



sin( 



t) = p^ sin 2/i -|- 2p{l — p) sin ^. 



(23) 



We assume a very large number of particles, so that the mean value of Eg. (p^) over different realizations of the spin 
lattice, is effectively obtained by considering only a single realization consisting of a very large number of particles. 
The enumerator is obtained simply as products of terms like in (|23|). The denominator of Eq.(16) involves the mean 
value of the square of a sum of products of cos- functions, and in the limit of large TV, it is dominated by terms which 
are products of four uncorrelated cos-functions, and we obtain 



l + (l-p(l 



As shown in Fig. 3 and observed already in 



COSfl))^ p^ SIT? ^ 



(r? sin 2fi + 2p(l — p) sin ^) 



(1 - cos^)) 

spin squeezing is obtained also in the case of a random filling. 



(24) 
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FIG. 3. Squeezing parameter as function of time in the Ising model with random coupling coefficients. Every coupling 
coeffecient Xi is a stochastic variable, attaining the value x with probability p and the value with probability (1 — p). The 
curves are obtained with the following values of p: 0.25 (solid curve), 0.5 (dashed), 0.75 (short-dashed), and 1 (dotted). 
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IV. CONCLUSION 



We have obtained general expressions for squeezing in different models of interacting spins. In the general model 
considered in section |ll| we showed that almost any kind of spin-spin interaction can be used to create squeezed states. 
As long as the coupling matrix in Eq. (||) is not proportional to the identity, spin squeezing can be produced. In the 
one-dimensional Ising model with arbitrary nearest neighbour coupling constants we obtained analytical results for 
all times, and we addressed the case of constant couplings and different kinds of departure from constant couplings 
of the spins. We found that the amount of spin squeezing was larger in a dimerized model with periodically varying 
coupling coefficients than in a homogeneous model. 

Among many possible generalizations of the present work we wish to mention long-range interactions and interac- 
tions in higher dimensions, and, e.g., the question whether extensions of the dimerized model to these cases are also 
superiour over homogeneous couplings. 

As mentioned in the introduction, spin squeezing is both a useful property in itself |^ and a signature of entanglement 
[^,0. It can be accomplished and detected in different physical systems without the need for experimental access 
to individual spins, and as such it can be used to test, e.g., spintronics proposals for quantum information processing. 
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